Magnetization plateaus and sublattice ordering in easy axis Kagome lattice 

antiferromagnets 
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We study kagome lattice antiferromagnets where the effects of easy axis single-ion anisotropy (D) 
dominates over the Heisenberg exchange J. For S > 3/2, virtual quantum fluctuations help lift the 
extensive classical degeneracy. We demonstrate the presence of a one-third magnetization plateau 
for a broad range of mag netic fields J 3 /D 2 < B < JS along the easy axis. The fully equilibriated 
system at low temperature on this plateau develops an unusual nematic order that breaks sublattice 
rotation symmetry but not translation symmetry — however, extremely slow dynamics associated 
with this ordering is expected to lead to glassy freezing of the system on intermediate time-scales. 

PACS numbers: 75.10.Jm,75.10.Dg 



Geometrically frustrated magnets, which are charac- 
terized by a large number of symmetry unrelated classi- 
cal ground states, display a wealth of new phenomena. 
In contrast to their unfrustrated counterparts, in which 
the low temperature ordering is largely determined by 
classical energetics even for small spin length S, quan- 
tum effects can play a crucial role in frustrated magnets, 
precisely because the classical energetics fails to pick a 
unique state. Of particular interest are those systems 
in which a broad cooperative paramagnetic regime [l| at 
intermediate temperatures gives way at low temperature 
to a variety of novel ordered and liquid phases Q aris- 
ing from the quantum fluctuations of spins. Even when 
magnetic order is selected, the resulting patterns are of- 
ten complex, and are promising candidates for realizing 
multiferroic properties Q. 

The Kagome lattice composed of corner sharing tri- 
angles is one of the the most frustrated lattice arrange- 
ments possible, and one that occurs commonly in na- 
ture. Ground states of the antiferromagnetic Heisenberg 
model on this lattice typically involve coplanar arrange- 
ment of spins, and are extensively degenerate in the limit 
of classical spins Despite much work on the role of 
quantum fluctuations in lifting this degeneracy, several 
fundamental questions regarding the ground state of low 
spin Kagome antiferromagnets (eg. S = 1/2, 1, 3/2 . . . ) 
remain to be decisively settled [2J. 

In some cases however, collinear spin arrangements are 
preferred, for example in the presence of strong single-ion 
anisotropy that leads to an easy axis. This occurs in the 
recently studied Nd 3 Ga 5 SiO 14 (NGS) compound 0, @, 
where Nd 3+ (total angular momentum J = 9/2) ions 
form a kagome antiferromagnet in which the common c 
axis is the easy axis at low temperature (below 33-ftT). 
The classical ground states of such collinear spins are 
also extensively degenerate. In contrast to the isotropic 
case, the effect of quantum fluctuations in selecting the 
low temperature state can be studied in a much more 



controlled fashion, and this is the focus of our work here. 

Our results are readily stated. First, the selection 
mechanism depends strongly on the spin length S. Be- 
low a critical spin length S c = 3/2, quantum transitions 
between different classical ground states dominate the 
degeneracy splitting. This 'kinetic energy' dominated se- 
lection applies only to S = 1 (5 = 1/2 does not admit 
a single ion anisotropy): There, a uniform spin nematic 
ground state was found in zero field, and a one-third 
magnetization (m = 1/3) plateau with ^/Sx v3 collinear 
order obtains in the presence of a magnetic field directed 
along the easy axis [7( . 

As we show here, the situation for all S > S c = 3/2 is 
entirely different. In these cases, virtual quantum tran- 
sitions dominate the energetics, leading to 'potential en- 
ergy' differences between different classical ground states. 
The m = 1/3 plateau, which is found to exist over a 
broad range of fields, displays an unusual nematic order 
at low temperature, i.e., it spontaneously breaks sublat- 
tice rotation symmetry but not translation symmetry (or 
spin rotation symmetry about the field direction). Fur- 
thermore, the system exhibits slow glassy dynamics in 
this state as a consequence of the free energy landscape 
induced by the 'potential energy'. 

In zero field, we obtain a simple characterization of the 
collinear states selected by virtual quantum transitions. 
The eventual low temperature behaviour in zero field de- 
pends on further entropic effects, and poses an interest- 
ing open question that will be addressed separately 0. 
These results have direct experimental relevance for the 
low temperature, finite field state of high spin Kagome 
magnets with easy axis anisotropy, such as NGS p, @] j 
and perhaps other collinear Kagome magnets as well. 

Model: Consider a layer of Kagome lattice antiferro- 
magnet, in the limit where the easy axis anisotropy dom- 
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inates over Heisenberg exchange, with a 'spin' S > 3/2: 



H = Si ■ - D £(S?) 2 -B^Sf, 



(1) 



(ij) 



where J > denotes the nearest neighbour antiferro- 
magnetic spin exchange interaction between the S > 3/2 
ions, D > is the single- ion anisotropy that picks out 
the common easy axis z, and B, the external magnetic 
field along this axis, has been scaled by the magnetic mo- 
ment gfJ-B- Due to the frustrated nature of the exchange 
J, the anisotropy effects can begin to dominate and pick 
collinear states for not- very-large D / J . In this collinear 
regime, we expect an analysis based on the smallness of 
J/D to give reliable results. 

With this in mind, we use J / D as the small parameter 
in a systematic perturbative approach that allows us to 
calculate the effective low-energy Hamiltonian and the 
resulting low-temperature phases. While our focus here 
remains the kagome case, we note parenthetically that 
our methods generalize readily to closely related models 
of triangular lattice magnets in zero magnetic field; these 
will be discussed separately [8|. 

Method: Our analysis proceeds by splitting the Hamil- 
tonian as H = H + XHi, where H = Jj2{i 3 ) S?S* — 
DYsiiSI? -BE, Sf and H x = ' E {lj) (S+Sf + h.c). 
Here Sf — Sf ± iSf, and A is introduced as a book- 
keeping device (A is set to one at the end of the calcu- 
lation). In each case discussed below, we begin by using 
standard degenerate perturbation theory Q in A to ob- 
tain the low-energy effective Hamiltonian that encodes 
the 'slow' dynamics induced by the Hi term within the 
ground state manifold of Hq. Terms in this expansion 
are naturally organized according to the power (n\) of A 
they carry, and by the number {rib) of different bonds on 
which constituents of H\ act. 

As it is possible to obtain a precise characterization 
of the pairs (n\,rib) that contribute to leading orders 
in the physical expansion parameter J/D, this proce- 
dure leads us directly to the physical low-ener gy e ffec- 
tive Hamiltonian at leading orders in J/D. [Ref [Uj em- 
ployed a closely related procedure for collinear states of 
the pyrochlore antiferromagnet in a field, but without 
anisotropy, and our calculation is better controlled due 
to the presence of anisotropy.] 

The kagome magnet at to — 1/3: We begin by consid- 
ering values of field such that the ground state manifold 
of Hq has magnetization to = 1/3 and is characterized 
by a 2 : 1 constraint that requires two spins in each tri- 
angle to be maximally polarized along the field, and one, 
minority, spin in each triangle to be maximally polarized 
anti-parallel to the field. While the Zeeman energy gap 
that drives the formation of this magnetization plateau is 
largest for B ~ 2 JS, we will argue later that the plateau 
extends down to a relatively small onset field that scales 
as B onset ~ J 3 /D 2 . 
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FIG. 1: (color online), (a) (Sublattice) rotation symmetry 
breaking on the m — 1/3 plateau: Dotted lines denote alter- 
nating arrangement of spins with average moment zero; up 
arrows correspond to S z ~ +S. c) First order transition to 
ordered state. 
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FIG. 2: (color online). Data for single-spin autocorrelation 
time t fit to the Vogel-Fulcher form (log(r)) -1 = (T-T/)/A; 
the intercept on the x axis yields Tf, while the slope gives 
A" 1 . 



A convenient way to represent the ensemble of to = 1/3 
states [ll[ is to encode the presence of a minority spin 
by placing a dimer on the corresponding bond of the 
underlying honeycomb lattice (see Fig [T] b). The low 
temperature physics on the plateau is then determined by 
the leading order effective Hamiltonian that acts within 
this dimer subspace. 

Because of the strong 2 : 1 constraint, the first term 
(for any S > 3/2) that breaks degeneracy of states in 
the dimer subspace is a diagonal (potential energy) term 
that occurs at order J 6 /D 5 ; this term arises from rc.& = 6 
n\ = 6 processes in the Van-Vleck expansion in which the 
six bonds on which Hi acts form a hexagonal loop. In 
contrast, the leading off-diagonal term, that corresponds 
to changing the state of a flippable hexagon (with alter- 
nating ±S values of S z around the hexagon) by reversing 
all spins on it, arises at order j6S-2^£)6S-3 ( w j^ n nfc _ 3 
n\ = 6S). 

A careful evaluation of these contributions for S = 3/2 
yields the effective Hamiltonian: 
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C1J 6 



c J 7 

|Ol)<Oi| - ^(\°3a)(o 3 b\ +h.c.) 



where O m (to = 0,1,2,3) denote hexagonal plaquettes 
with to dimers and O3A, 035 represent the two flip- 
pable dimer configurations with three dimers on a hexag- 
onal plaquette. The values of the coefficients are Ci = 
w 0.1335 and c 2 = ^ « 3.3 x 10" 3 . Thus, the 
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off-diagonal term of 7i m — 1/3 is negligible compared to the 
diagonal part for S — 3/2. From the foregoing analysis, it 
is also clear that its magnitude decreases with S, and we 
have therefore not calculated the off-diagonal coefficient 
for S > 3/2. The potential energy V(S) = c t (S) J 6 /6D 5 
has however been calculated for general S, and the result 
is V(S) = (2S*) 6 J 6 /1024(25- 1) 5 L> 5 . 

As an independent check on this result, we use the 
semiclassical large-S* expansion procedure of Ref [13] and 
directly calculate the semiclassical effective Hamiltonian 
in this large-S* limit: The leading 0(J 6 /D 5 ) term ob- 
tained by expanding the semiclassical result in powers of 
J/D agrees precisely with the large-S* limit of the per- 
turbative result obtained above for arbitrary S. 

The low temperature physics is thus well-described by 
a classical dimer model with weights associated with the 
potential energy term in the effective Hamiltonian. Fur- 
thermore, this classical potential energy is minimized by 
any configuration with no hexagon having precisely one 
dimer on it. As there are a large number of such config- 
urations, a more detailed analysis is needed to elucidate 
the nature of the low temperature state. To this end, 
we have employed a generalization Q of the procedure 
of Refs [l3l . to efficiently simulate an interacting clas- 
sical dimer model with this potential energy term — our 
algorithm employs non-local loop updates but preserves 
detailed balance in order to generate the correct equilib- 
rium Gibbs distribution at temperature T. 

As the temperature is lowered to below T c « 0.23V (S) 
(obtained by extrapolating the finite L data in FigfT]), 
we find that the system undergoes a transition to a state 
with (sublattice) rotation symmetry breaking as shown 
in Fig Q] a: In this simplest schematic of the ordered 
state, one spontaneously chosen sub-lattice of spins ac- 
quires the maximal polarization +S along the field. In 
order to satisfy the strong 2 : 1 constraint on each trian- 
gle, the spin moments on the other two sub-lattice sites 
then alternate +S, — S . . . in one of the two possible al- 
ternating arrangements along a stripe. [A very similar 
ordering was suggested earlier for an isotropic kagome 
magnet in the semiclassical limit 15[ .] This ordering can 
be conveniently characterized using the sublattice order 
parameter $ = ^2 p m p e 2p7 ™/ 3 , where m p denotes the sub- 
lattice magnetization of the p th sublattice (Fig[T]a). The 
two peak structure in the histogram of |$| 2 at T c provides 
evidence for the first-order nature of the transition. 

As each stripe can be in one of two possible alternat- 
ing states, its internal state can be represented by an 
(Ising) pseudo-spin variable a. Are these a ordered in 
any manner, or do they fluctuate over time? The absence 
of any non-zero wavevector bragg peaks in the numeri- 
cally measured static structure factor of physical spins 
rules out any q 7^ order. Furthermore, the low temper- 
ature value of |$| 2 is, within error bars, exactly what one 
would expect if each stripe fluctuated between its two al- 
lowed alternating states (ruling out q = order for the 



pseudospins). Additional confirmation also comes from 
the statistics of different types of hexagons with precisely 
two dimers on them (not shown). 

This throws up a point of general interest: Although 
the fully equilibriated system only breaks rotation sym- 
metry by forming disordered stripes, this equilibriation is 
achieved in our numerics because the algorithm incorpo- 
rates non-local loop updates that can flip a macroscopic 
number of spins in one move. In the experimental sys- 
tem, the dynamics is of course purely local. Such local 
spin flips cost significant potential energy and the system 
needs to change the internal state of an entire stripe to 
avoid this potential energy penalty. Systems with very 
similar potential energy landscapes have been the sub- 
ject of earlier studies which demonstrate that the time- 
scale for changing the internal state of a stripe diverges 
rapidly with system size if the dynamics is local [l6| . It is 
thus clear that the low temperature phase displays glassy 
freezing of the stripes. 

We have also numerically studied the behaviour of 
the system for T c < T with purely local two-spin ex- 
change dynamics satisfying detailed balance, and moni- 
tored the temperature dependence of the single-spin au- 
tocorrelation function for a range of moderately large val- 
ues for the ratio JS 2 /V(S). In these simulations, the 
magnetic field B is fixed to its nominally optimal value 
B = 2JS which places the system close to the center of 
the m = 1/3 plateau; however, we emphasize that the 
ratio JS 2 /V(S) is kept finite as we wish to explore the 
higher temperature dynamics, and configurations outside 
the 'dimer subspace' are allowed but exponentially un- 
likely (as opposed to forbidden). We find that the single 
spin autocorrelation time r increases very rapidly as we 
lower the temperature. Indeed, r(T) can be fit well by 
an activated functional form of the Vogel-Fulcher type 
t(T) = exp(A/(T-Tf(L)) (Fig©, thus extending fur- 
ther the analogy to other models of glass-formers [l(| . In 
our fits, the freezing temperature Tf(L) drifts somewhat 
with linear size L, but its extrapolated L — > 00 value is 
within 10% of the equilibrium T c obtained earlier, while 
the barrier energy scale A shows no L dependence. 

We thus conclude that while the fully equilibriated sys- 
tem only breaks sublattice rotation symmetry but not 
lattice translation symmetry, slow glassy dynamics that 
sets in as the temperature is lowered through T c forces 
the system into a glassy state in which the stripe pseu- 
dospins also freeze in a random pattern, thereby break- 
ing lattice translation symmetry in a random manner. 
What would be a good experimental signature of this be- 
haviour? Clearly, it is not appropriate to focus on elastic 
Bragg peaks within the first Brillouin zone, as neither 
the fully equilibriated sublattice ordered state nor the 
metastable glassy state lead to such a Bragg peak. How- 
ever, we note that the sublattice order parameter that 
indicates the breaking of rotation symmetry can be re- 
constructed by taking suitable linear combinations of the 
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FIG. 3: (color online). The three types of unfrustrated bonds; 
presence (absence) of red dots represents S z = +S (S z = — S) 



measured spin structure factor close to wavevectors with 
components (0,0), (0, 2n/a), (2n/a,0), and (27r/a, 2-k/cl) 
along To and Ti (FigQ]a), and this provides a possi- 
ble experimental probe of sublattice symmetry breaking, 
that is independent of more subtle questions regarding 
the freezing of the stripe pseudospins. 

The kagome magnet in zero field: To understand how 
the onset field of the m = 1/3 plateau scales with J/D 
and S, we need to first characterize the collinear states se- 
lected by analogous potential energy effects in zero field. 
We begin by noting that at large D and B = 0, the 
ground states of Ho are obtained by requiring all spins 
to have polarization ±S along the z axis, and allowing 
only one frustrated bond (pair of nearest neighbour par- 
allel spins) per triangle. 

This macroscopic degeneracy is broken by the effect 
of real and virtual quantum transitions. A perturbative 
analysis in J/D again allows us to derive a leading or- 
der effective Hamiltonian that acts within this degenerate 
subspace and encodes these effects: 

Hb=q = -W^(2|3B)(3B| + |2B)(2B| + 0|1B)(1B|) 

B 

+ t^2{\W){W\+h.c.) (2) 



m = 1/3 plateau. Thus, the onset field will scale as 
Bonset ~ J 3 /D 2 . Since this onset field does not scale 
with S (and can be quite small for even moderate values 
of D/J), we conclude that the predicted magnetization 
plateau state is likely to fall well within the field regime 
accessible to experiment even for large S. 

Discussion: We have thus predicted an unusual sub- 
lattice ordered m = 1/3 magnetization plateau state 
with slow glassy dynamics at low temperature in pure 
S > 3/2 Kagome antiferromagnets with strong easy axis 
anisotropy. We have also provided a simple characteri- 
zation of the collinear states that are selected in zero ex- 
ternal field by an effective potential energy generated by 
virtual quantum fluctuations. Our results are expected 
to provide an excellent starting point for understanding 
the very low temperature physics of the kagome anti- 
ferromagnet NGS, although more work may be needed, 
particularly in the zero field case, to understand the 
effects of spatial distortion and sub-dominant in-plane 
anisotropies [H . We hope our results provide impetus for 
experiments on other easy axis kagome antiferromagnets, 
especially members of the large family of materials [B[ 
that have the Ca3Ga2Ge40i4 crystal structure of NGS. 

We would like to acknowledge useful discussions and 
correspondence with L. Balents, B. Canals, R. Cava, 
D. Dhar, R. Moessner, T. Senthil and F. Wang, com- 
putational resources of TIFR, and support from a Hell- 
mann Fellowship (AV), LBNL DOE-504108 (AV), and 
DST SR/S2/RJN-25/2006 (KD). 



where the sum extends over unfrustrated bonds B 

27J 3 ± _ 9J 3 



(ilD- 



For 

, 2] yz, \MB) represent 



S = 3/2 we find W 
unfrustrated bonds with local environment of type M 
(Fig (3j), and |3B) and |3B ) in the second (off-diagonal) 
term are related to each other by a spin-exchange be- 
tween the antiparallel spins connected by the flippable 
{i.e. type-3) bond B. More generally, for S > 3/2, the t 
term is negligible while W(S) = 2D^(2S-i)' 2 • 

Constraints provided by the kagome geometry allow 
us to prove that the potential energy term W is mini- 
mized by a class of collinear configurations in which no 
spin is the minority spin of both triangles to which it 
belongs 0] , and demonstrate that this class of configura- 
tions has macroscopic entropy 0. 

More importantly for our purposes here, this zero field 
ensemble is stable to small magnetic fields since the defin- 
ing constraint on minority spins does not fix the average 
magnetization. When the field is increased further be- 
yond B ~ J 3 /D 2 , the 0(BS) Zeeman energy gain of the 
rn = 1/3 ensemble (defined by the 2 : 1 constraint) will 
begin to dominate over the 0(J 3 S/D 2 ) potential energy 
gain of the zero field ensemble defined by the constraint 
on minority spins, and this will trigger the onset of the 
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